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ABSTRACT 


TIwce methods for evaluating the accuracy of hydrograrfhic 
positioning data are presented. One method consists of 
classifying each position in a survey based on the radius cf 
the 90 percent confidence circle. The second method 
involves classificaticn of positions based on tne farameters 
of the 90 percent confidence ellipse. Both methods are 
based on gecmetric and statistical relationships Lretween 
intersecting lines of position. 

Range-Lrange, azZinuth-azimuth, and range-azimuth [csi- 
tioning data are classified using both criteria. For 
noncritical positions, the confidence circle method is found 
to be preferable due to its ease of interpretation. For 
fositicns of significant features, such as underwater 
hazards, the confidence ellipse provides a more useful 
representation of the shape and orientation of the true 
error distribution. 

The concept of presurvey fositioning design is alsc 
presented. With the aid of computer graphic displays, the 
hydrograrher can predict the accuracy of offshore posi- 
tioning data prior to data acquisition. By analyzing accu- 
racy lobes generated about shore stations, a survey can be 


designed to meet given specifications. 
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I. INTRODUCTION 


A. FPACKCRCUND 


A hydrographic record can te viewed as the resultant cf 
two independent measurements made at a discrete point cver a 
body of water. These measurements involve the determination 
of a vessel's positicn at a given time as well as the depth 
of water at that position. Of interest to the hydrograrher 
and tc the user of hydrographic data is the accuracy of tte 
position determinaticrs. Fundamental to the determination 
cf positional accuracy is the identification of the sources 
of errors in positicn measurements and the ultimate treat- 
ment cf these errors. 

A hydrcgraphic position can be determined by a number of 
methods all involving geometric relationships betweer known 
points and the vessel's unkncewn location. The known foints 
may be fixed stations on shore, whose coordinates have Leen 
determined by geodetic survey methods, or they may be 
rapidly roving sateliites whose coordinates in time and 
space can ke defined very precisely. A hydrograrhic pcsi- 
tion is established Ey the intersection of two or more lines 
of position (LOP's) which are generated by the geometric 
relaticnshirs between the fixed points and the vessel's 
unknown location. The resultant accuracy of the vessel's 
position is therefore, in part, a function of the errors 
associated with the irtersecting LOP's. 

Séveral neasures cf accuracy can be used to evaluate the 
quality of a hydrographic position. Predictability, or 
absolute accuracy, is the measure of accuracy with which the 
positioning system can define the location of the same foint 


in terms of geographic coordinates. Repeatability, cr 


relative accuracy, is a measure with Whten a postertornieg 
system permits a user to return to a specific point on the 
earth's surface in terms of the LOP's generated [Ly tne 
system (Ref. 1, p. 14j. With the elimination of all systcn— 
atic cr kias errors, the terms repeatability and predicta- 
rility become identical. Hydregraphic surveyors usually 
work toward this condition, although it is not always 
achievakle. 

Heinzen [Ref. 2] and Burt [Ref. 3] have presented 
several techniygues fcr quantifying the repeatable accuracy 
for offshore positions. These technigues have roots in tte 
statistical treatment of randcm error. Aithough the methcds 
have keen well documented, no single criterion to classify 
the accuracy of a hydrographic position has been agreed ufon 
ky the internaticnal hydrographic community. 

Freceding the development of automation in hydrograrhic 
data acquisition and frocessing, the task of calculatirg an 
accuracy figure to attach to each position ina hydrographic 
survey was unthinkable. To ensure overall accuracy in a 
survey, certain generalizations were developed to act as 
guidelines. For example, the U.S. Coast and Geodetic Survey 
Hydrographic Manual [Ekef. 8, p.j 217] states the following 
concerning the strength of a three-point fix: 

The £1x is SE ond when the sum of the two uangiles vis 
egual to or greater than 180° and neither angle is less 
than 30°. The nearer the angles equal each other the 
Strenger will be thesrix 
Generalizations of this type provided useful qualitative 
guidance for assuring a degree cf positional accuracy and 
many are still in existence today. 

With the aid of ccmputers, the hydrographer now has the 
capacity to evaluate the accuracy of positioning data fcr an 
entire survey. An accuracy figure can be computed for each 


Fositicn in a survey and stored in a data base along with 
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cther survey information. This figure may provide useful 
information for users of the data, as well as a yardstick 
for the hydrograrher to evaluate the guality of the work. 
Furthermcre, a presurvey accuracy analysis enables a survey 


to be designed to meet desired specifications. 


EF. ACCURACY STANDARDS FOR HYDBOGRAPHIC POSITIONING 


In 1982, the Intérnational Hydrographic Organizaticn 
(THO) puklished new recommendations for error standards 
concerning the accuracy of hydrographic positions. These 
standards [Ref. 5] are: 

The position of sourdings, dangers and ali other signif- 
icant features should be determined with an accuracy 
such that any Frobakle error, measured relative to Shore 
contrcl, shall seldcm exceed twice the minimum flottarle 
error at the scale cf the Sune (normally 1.0 mM on 
paper). it is most desireabie that whenever positions 
are determined by the intersection of lines of fpositicn, 
three such lines be used. The angle between any fair 
should not be less than 30°. 

Mest statisticiars define the term “probable errer" as 
that e€xricr occurring at the 50 percent probarility level. 
However, the author cf the IHO Standards, Commodore A. #. 
Cooper RAN (Ret.) has stated that the term “probable error" 
waS intended to nave no statistical significance. Munscn 
interpreted the words "shall seldom exceed" to mean 10 
percent of the time [Fef. 6]. Using this interpretaticr, 
the first sentence of the specification might be written: 

The bes eeton of soundings, dangers and ail other signif- 

icant features should be determined with an accuracy 

such that any error in_ position measured relative tc 

Shere control will fall within a circle with radius of 

the wirimum plottable error at tae scale of the survey 

(normally 1.0 mm. on paper), with 90 percent confidence. 
The specification in this form could be evaluated gquantita- 
tively. The critericn for defining accuracy in terms of a 
fixed probability is common in the field of surveying. For 


example, the standards of accuracy developed for geodetic 
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contrcl surveys have their origin in probability thecry. 
Procedures for obtaining first-crder geodetic positions 
require sixteen repeated theodolite observations of each 
direction. Lower order positions reguire fewer numbers of 
cbservations. Given the precision of one observation cf 
each direction, it can be demonstrated that increasing the 
rumber of observations coincides with increasing the frcb- 
ability cf the direction falling within specified limits. 

Regarding accuracy determinations, there are several 
Froblems unigue to hydrographic surveying. Whereas stan- 
dards for other types of surveys rely on multiple observa- 
tions of the same guantity, the accuracy of a hydrograrhic 
position must be evaluated in terms of a single observation 
(which may be the intersection of two or more LOP'‘s). 
Diverse methods for oktaining a hydrographic position exist 
and these methods must all be evaluated using the same 
criterion. Also, there is a broad spectrum of eguipment 
used in hydrographic fositioning and in many cases the 


precision of this equipment is not well defined. 


C. CEJECTIVES 


A need exists to give guantitative meaning to the accu- 
Lacy specifications set forth by the IHO. One of the okjec- 
tives of this thesis is to demcnstrate that defining the 
specifications in terms of the fixed 90 percent confidence 
levei 1s a valid interpretation. By defining what the spec- 
ificatiors imply, procedures can be developed to meet the 
standards. 

A second objective of this thesis is to apply tke thecry 
of @C€L©rcrs, associated with hydreqca phicwpocie cn 11g toma 
data set. This analysis involves classifying positicning 
data acquired in a survey based on the radii of circles of 
equivalent probability. It will be demonstrated senate es 


method of classification is a useful index for guantifying 
the accuracy of fositions. The computed radii of the 90¢ 
percent confidence circles can serve as an accuracy figure 
that can ke attached to each position in a survey and stored 
in a cata base. 

The third objective of this thesis is to demconstrate 
that a presurvey analysis can be used in designing pcsi- 
tional accuracy to meet specifications. The existing 
general guidelines for planning can be tetter defined. For 
example, in planning a survey hydrographers usually lay out 
circles which delimit the 30° and 150° roundaries that 
define the minimum and maximum allowable intersection angles 
ketween two LOP's. As a means to meet accuracy require- 
ments, it can be shown that these limits should vary kEased 
on the scale of the survey and the precision of the fosi- 


tioning eguipment. 
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II. NATURE OF THE PROBLEM 
The development cf an accuracy figure for offshcre fosi- 
tions is inherently tied to the geometry of the fositionirg 
methcd and the errors which are associated with the fosi- 
tioning eguipment that is used. This chapter will discuss 
the gecmetric and statistical elements involved in deéeter- 
fining an offshore position and presents several methods for 


guantifying repeatable accuracy. 


Aw HYDRCGRAPHIC POSITIONING GEOMETRIES 


An offshore £1ix can be determined by the intersecticn cf 
two or more LOP's. These LOP's may be generated by elec- 
tronic or visual means. Working toward the development of 
an accuracy index, it will be necessary to compute the ancle 
of intersection of tke LOP'S associated with different fosi- 
tioning geometries. The following sections discuss the 
geometry cf conventicral offshore positioning methods and 
ways to compute the argles of intersection. This thesis 
wiil not address the geometry involved in a three-point 


eextant f1x. 
1. Range-Range 


Fstablishing an offshore fix by range-range gecnmetry 
involves measuring distances electronically from fixed fosi- 
tions on shore to the vessel's unknown location. Ranges can 
ke determined by measuring the elapsed time between trans- 
misSion and receipt of a radio pulse or by comparing the 
Fhase of the transmitted wave with the phase of the received 
wave [Ref. 2]. In each case, transmitters are set on 
stations on shore whose coordinates are determined by 


Frecise land survey methods. 
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An electronic pcsitionirg system may te active cr 
passive. In an active system, a transmitter from the survey 
launch keys the transmission of ranges from the shore 
Staticn. In turn, the signals generated from the shcre 
staticns (slaves) are then received by the launch. An 
active system is limited to a finite number of users, 
usually not more than about four. The number of users cf a 
passive system is unlimited as the survey launch reguires 
only a receiver which is constantly listening for signals 
which are EFeing transtitted from shore. 

Short-range, cr line-of-sight, positioning systems 
are used for nearshore hydrographic surveys. These systems 
cperate in the microwave region of the electrcemagnetic spec- 
trum (3 to 10 GHz). A distance is determined by observing 
the time needed for a fulse to travel from a master tran- 
sponder located aboard the survey vessel to a remote tran- 
Sponder cn shore and back to the master transponder. 

Knowing the average velocity of the electromagnetic pulse, 
the distance D is theLz 


| 2. ico 


where c 1S the group velocity of the wave packet and t is 
the two-way travel tine. Shcrt-range systems which are in 
wide use today are Racal Decca's "Trisponder" and Motorcla's 
"Mini-Ranger." These systems have direct range readout and 
are readily interfaced into a navigational computer anda 
data acquisition system. Both systems are active and user 
iimited. 

Medium-range positioning systems operate in the 1- 
to 5-MHz frequency range of the electromagnetic spectrunr. A 
distance is determined by measuring the phase relaticnship 
ketween transmitted and received waves. These systems are 


usually referred tc as continuous wave systems and the 
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problem cf lane ambiguity must be addressed. Ranges are 


expressed in full and partiai lane counts ™where a anc jaca 


wis 


x 
i 
n> 


beZiaiey) 


where X is the wavelength of the transmitting frequency, f£, 


and given by 


(2.3) 


>» 
a 
+] 


Medium-ransge systems ccmmonly in use today are Cubic 
Western's "ARGO," Hasting Raydist's "“Raydist," and Odon 
Cfishere'’s "Hy drotrack.* 

The angle of intersection associated with a range- 
range poSition is computed frcm a simpie trigonometric rela- 
tionship. The vessel's position P (Fig. 2.1) 1s determined 
by the intersection cf the ranges from the left and right 
shore stations, R1 and R2 respectively. 38 is the base line 
distance computed between the two known shore stations. 
Since the range circles from the shore stations intersect at 
two points, it is necessary for the flotter to reccgnize 
which side of the base line the vessel iS on in order to 
eliminate the ambiguity. The angle of intersection cf the 
two LOP's (8B) is given by the law of cosines 


B?-R17-R2? ) 
eR ho 


B = 180° - Arc cos (2.4) 


In gualitative terms, the fix is strongest when 8 approaches 


90°. Most hydrographic specifications limit the angle cf 


intersection from a pinimum of 309 to a maximum of 1£09. 
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Figure 2.1 Gecmetry of a Range-Range Position 


Hydrographic positioning by hyperbolic-hypertolic 
geometry utilizes the intersection of two hyperbolas each 
generated about a pair of Shore ccntrol stations. A hypfer- 
bola is the locus of points in which the difference cf 
distance from two fixed points is always constant. A three- 
staticn hyperbolic net is the most commonly used hypertolic 
mode fcr offshore survey (Fig. 2.2). One family of hyrer- 
rolas (Red) are generated about a master station, W, and a 
Slave, R; while a seccnd family of hyperbolas (Green) are 
generated with respect to the master and a second slave, G. 
bere tne finst fLamiiy cf hyperbolas, the control points ™M and 
R act as the foci, while points M and G act as the fceci fer 


the seccnd family. 
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Hyperbolic location methods can be divided into two 
groups based on the electronic principlesmused sto gdetrnewiae 
distance differences [Ref. 7, pp. 87]. Loran is an example 
of a pulse system in which the differences in times cf 
artival of pulses transmitted Fy the master-slave combina- 
tions are translated into distance differences. The resul- 
tant position has no lane ambiguity and is easily resolved. 
The seccrd method of hyperbolic positioning involves meas- 
uring a phase differerce from two master-slave combinations 
at the vessel's position. The phase difference translates 
into a fractional lane count which in itself provides an 
ambiguous position. This ambiguity is resolved Ey using a 
whole-lane counter which is initialized at a known geografh- 
ical point. In hyperbolic positioning, the. Sum ots eera 
passive mode and the system can be used by many vessels. 

The angle of intersection between the two hyrertoias 


can be ccmputed by first defining the following quantities: 


S. is the length of red base line, 

3 is the length of green tase line, 

RS is the distance tetween master and vessel's pcsiticn 3D, 
R. is the distance from red slave to point P, 

R, is the distance from green slave to point P, 

OL is the angle between lines PM and PR, and 

o is the angle between lines PM and DG. 


The spacing Fetween lanes increases with distance 
from the master-slave pair. The lane widths along the Lase 


line are 


d and w! X 
g = Peas 
5 = ( ) 


Then the lane widths at any pcint P are 
d h l 


= l and w= 
"te aatasey ?  sintafay) (2-6) 


1é 


where tke term 1/sin(a/2) is called the lane expansion 
factor. The angle of intersection ,8, between the twe 
hypertolas is then given by 


ae ee (27) 
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Figure 2.2 Gecmetry of a Hyperbolic-Hyperkolic Position 
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- Range~Azinuth 


This positioning geometry is used for nearshore, 
line-cf-sight surveys. One LOF is generated by an elec- 
tronic range originating from a transmitter located on a 


shore ccntrol station. A microwave system is commonly used 


i 


in this arrangement Lut systems employing a laser can alsc 
be used for short-range work. Another LOP is generated by 
fixing an azimuth frcp a shore control station to the 
vessel: A second control station 15 90secdmror agen Gras 
azimuth Ly the observer. Azimuth determinations can be made 
after cbserving directions with a theodolite as an observer 
tracks the moving vessel. 

There are twe ways to determine a range-aZinuth 
positicn. The most ccmmon way is to have the theodolite and 
the transmitter occupy the same shore control station. 
Hence, the angle of irtersection, g, of the LOP's is always 
S09. This arrangement 1s commonly used by the National 
Ccean Service {NOS) fcr large-scale nearshore surveys. 

The other way is to have the theodolite and the 
transmitter occupy two different control points. Then the 
geometry iS Similar to that of the range-range position. 

The angle of intersection, pg, 1S computed by trigoncmetric 
relaticnshifs among the azimuth of a line between the skcre 
Staticns, the observed direction to the vessel, and the 


measured range to the vessel. 


YW. Azimuth-Azinuth 


> ee ee ee ee ee ee ee ee 


AZimuth-azimuth positioning geometry is used for 
nearshore high-accuracy surveying. Theodolites are set over 
two ccntrcel stations cn shore. The vessel is sighted on 
Simultaneously by the two theodolite observers, generating 
two visual LOP's whose intersection define the vessel's 
location. Initial azimuths are fixed by sighting cn centrol 
stations which are visible to the observers. 

The angle of intersection for an azimuth-aziruthk 
position is dependent on the geometric relationships ktetween 
the cccuried stations, the initial stations, and vessel's 
position (Fig. 2.3). Assuming that theodolite observers 


20 


Scemmeyreccat2ons | and 2, and anitial on stations 3 and 4, 
respectively, the observer at station 1 measures angie y, and 
the observer at staticn 2 measures y to the vessel. The 
angle of intersection, 8, is then poapured Byefirse deter- 
Mining the forward azimuths, measured clockwise from the 
Semen, <r¢cm stations 1 to 2 (ode 1 tom (a). 2G sl 
Cele and 2 to 4 oe The interior angles, 2 en oes of 
triangle 12P are 


oa lo ae - a, ae) 
and 
9 = |a + - - 
2 | 24 Me Oo (ez 2.9) 


so the angle of intersection, B , at the vessel's location 


is 


= 0 =~ 
8 = 180 (8 + 8) (2. 10) 


Eee CRASSES OF ERRORS 


All hydrograrhic positioning measurements are subject to 
error. The following sections discuss categories of errors 


and methods used to treat these errors. 


1. Blunders 





Plunders are aross mistakes which are generally due 
to the carelessness of the observer. Blunders can vary in 
Magnitude, ranging from large errors which are easily 
detected, to small Errors which may be barely distinguished. 
They can be detected by making repeated observations cr by 
carefully checking the data in the processing phase. 
Blunders occur in various forms and most can te avoided by 


carefully planning tte data acguisition process. 
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Figure 2.3 Geonetry of an AzZimuth-Azinuth Position 


Consider the following as an example of a blunder 


Pe ee a= 





associated with range-range geometry. An offshore fesiticn 


is to ke determined Fy the intersection of two electronic 
LOP's generated from transmitters located on known shore 
staticns. The vessel is working west of a shoreline that 
runs generally in a north-south direction. As the hydrc- 
grapher faces the stations from sea, the southern shore 
station is mistakenly identified as left and the northern 
shore station as right. The resultant offshore positicn 
will plot to the east of the base line. This blunder is 
readily detected and can be easily remedied. 

Not all types of blunders are so easily detected. 


Suppose an offshore pcsition is to be determined ry a 


rangé-azimuth fix. A range and an azimuth are generated 
from a kncwn control station to the vessel's position. A 
second control staticn is used to fix the initial azimuth; a 
third shore control station is located 10 meters from the 
initial station and its coordinates are mistakenly used for 
the initial station in plotting. The resultant hydrographic 
position is in error, but this error will not be easily 
distinguished. 

Although most blunders have their origin ir human 
carelessness, some can be attributed to equipment malfunc- 
tion. Fcr example, microwave systems which generate LOF's 
are known to become unsteady under certain conditior#rs. 
Spuricus range readings resulting from signal reflections 
can Le recerded as true positioning data. In this case, the 
blunder may or may not be easily detected. 

In automated data acquisition systems, software has 
keen developed te detect the occurrence of anomalous range 
readincs. By inputting a course and speed of a vessel trav- 
e€ling along a line, tke computer can determine if the 
recorded position is valid based on the principle of dead 
reckoning. If the recorded position is found to be invaiid 
the hydrographer will be immediately alerted to the situ- 
ation and can take action to remedy the problem. In non- 
automated systems the principle of dead reckoning is applied 
manually. Given the course and speed of the vessel, the 
validity of the position can be checked with spacing 
dividers. This involves checking the spacing between fixes 
recorded before and after the position in question. 

Pefore any tyfe of error analysis is to Le fperfcrned 
con the hydrographic fesitioning data, it is essential that 
all blunéers be identified and properly treated. In 
general, careful planring coupled with thorough checking 


will givdinize the occurrence of blunders. 
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2. Systematic Errors 


Systematic errers occur with the same sign, usually 
of similar magnitude, and can te expressed in terms of a 
mathematical model. Systematic errors follow a defined 
pattern and occur in a number cf consecutive related ckser- 
vations. Repetition cf measurements does nothing tc miri- 
mize their effect. In the case of hydrographic positioninxg, 
systematic errors are identified and modeled by calikraticn 
of the measuring instrument against a known standard. the 
following is a brief discussicn concerning systematic errors 
and their treatment in relation to hydrographic positionirg 


equipmert. 
a. Theodolites 


In néarshere surveys the theodolite is used 
Frimarily for range-azimuth and azimuth-azimuth positicning. 
Systematic errors asscciated with the theodolite can he 
classified into two groups: those associated with tke phys- 
ical design of the instrument and those involving the geon- 
etry of the positioning scheme. Some sources of systematic 
errors [Ref. 8] associated with the physical characteristics 
cf a thecdolite are: 
i. The horizontal circle may be eccentric. 
41. Graduations on the horizontal circle may not Le 
uniform. 
1ii. The horizontal axis of the telescope (about which 
1t rotates) may not be perpendicular to the 
vertical axis of the instrument. 
iv. The longitudinal axis of the telescope may net te 
normal to the horizontal axis. 
v.- The telescope axis and the axis of the leveling 


Eukble may nct be parallel. 
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These errors are usually small in magnitude and can be elin- 
inated by proper adjustment of the instrument by either tke 
Manufacturer or a qualified technician. 

The field hydrographer has ultimate control over 
the gecmetric systematic errors associated with a theodc- 
lite. In range-azimuth positioning the theodolite and 
transwitter may cccury the same horizontal control station. 
If the theodolite is not set directly over the staticn a 
resultant systematic error will occur in all measurements. 
it can be shown that these errors are non-linear but do 
follow a mathematical relationship. Likewise, if the trans- 
mitter is not located directly over the station, a similar 
type cf kias occurs. Depending on the eccentricity of the 
theodclite, the vessel's range from the theodolite, and the 
scale of the survey--these errors can seriously affect the 
absoiute accuracy of the offshcre positions. 

in a Similar fashion,,it is also,imperative tc 
positicn the target directly over the horizontal ccentrol 
staticn used aS an initial. Failure to do this wiil result 
in an error which will be propagated to offshore positicns. 

Many Situations arise in the field where it is 
advantageous to set a transmitter and theodolite over a 
ele Nemuzontameeontrol station. Frequently it is 
feasikle to construct a platform to accommodate both instru- 
ments; in a case where it is not, the position of an eccen- 
tric hcerizontal contrcl station near the original staticn 
should Fe determined and that station used for the lccaticn 
of one of the instruments. The theodolite and the trans- 
mitter then occupy txc known stations and the geometric 


source of systematic error is eliminated. 
b. Electronic Ranging Systems 


The systemfatic errors associated with electronic 


positioning systems are complex in nature and functicns of 


Many variables. Munscn {Ref. 9, p- 4] addresses several 
problems associated with short-range systems used ir hydro- 
graphic surveys. The most ccmmon problems with short-range 
systems are variation in range and calibration drift with 
time. Variations in internal equipment time delays in the 
transmitter, the transponder, or the receiver can induce 
errors in measured rarges. For pulse systems such varia- 
tions can occur due to temperature dependence of components 
and fluctuations in signal strength at the transponder. 
Multipath effects are also a problem. Under some circun- 
stances a reflected wave and the directly transmitted wave 
arrive with a phase difference of 1809. Cancellation cr 
fading of the directly transmitted Signal can result. 

NOS conducts base line calibrations of short- 
range positioning systems periodically during the course of 
a Survey to minimize cr eliminate systematic errcr. In this 
process, a transmitter and receiver are each placed over 
contrcel stations on shore and the measured range is ccmpared 
to the true range. Ir this way the systematic error is 
eliminated Ly zeroing the instrument or by applying a 
constant correction tc raw data. System checks are 
performed daily to assure there is no drift from the orig- 
inal calibration. A check can be accomplished by comparing 
a positicn defined by the ranging system to a known fixed- 
Foint position, to a sextant fix position, or arn intersec- 
tLOnmeeo sition, 

Munson [{Ref.j{ 9, p.- 5] also discusses sources of 
Systematic errors associated with medium-range systems. the 
most significant systematic errors occur as a function cf 
position due to varying propagation velocity. The mediun- 
range electronic signal propagation velocity depends on the 
surface conductivity and transmission path (over water, over 
dand, or over different types of land). Because of this 


dependence, systematic errors as a function of position 
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cccur at different effective phase velocities. Knowing the 
propagation velocity to use, or the phase correction tc nake 
as a function of rancse, is a problem. Sky wave and storm 
interference also pose problems. At extreme ranges of cper- 
ation, sky wave interference can affect the more predictable 
ground wave, especially during nighttime operations. lLane 
ambiguities are also a problem. Most systems are inherently 
ambiguous and must be zero set and continually mcnitcred for 
dane jumps or loss of signal which results in the loss cf 
eae Coun ti 

NOS uses several technigues to determine the 
systematic error asscciated with medium-range positionirg 
systems. These technigues involve determining a whole and 
partial liane count for phase ccmparison systems. Two of the 
more widely used technigues are comparison of three-foint 
sextant £f1x positions to positions determined by the elec- 
tronic ranging system and calikration of the electronic 
system at a fixed feint. In both technigues the whole lane 
counts are fixed by the calibration; correctors to the 
partial lane count are determined and applied to the raw 


ranging data. 


Random errors are chance errors, unpredictable in 
Magnitude or sign, and are governed by the laws of frob- 
ability [Ref. 10, p. 1206]. They are errors which remain 
after blunders and systematic errors have been removed. 
Randcr errors result from accidental and unknown ccmbina- 
tions of causes and are beyond the control of the observer. 
CheCiweimieretem—lz, 6 b= 2 estates they are Characterized by; 

ij. Variation in sign; positive errcrs occur with 

egual frequency as negative ones. 

ii. Small errors cccur more frequently than large errors. 


iii. Extremely larce errors rarely occur. 


oa 


Random errors are unigue to specific types of fesi- 
tioning eguipment and vary in pagnitude depending on the 
precision of the instruments that are used. The following 


secticn cutlines statistical methods for their treatment. 


Cw. TREATMENT OF RANICA ERRORS. 


1. Cne-Dimensional Errors 








Certain basic statistical quantities must first be 
defined in the analysis of random errors. Consider a vessel 
moored securely to a fixed offshore platform. A numter of 
ranges, n, from a microwave transmitter located on a shore 
control station are recorded. The mean of these okserva- 
tions is 

as ; <i 
{zl (Zonta) 
where x, represents an individual observation. The standard 


error, Ss, of the observations is then 


n 
2 - 
> ene om (x, a uy) (22°42) 


where the guantity (2-H ) is referred to as the residual, 
cr true €rror, Vis ofa particular observation. As n gets 
very large, the factor 1/n can be substituted for 1/(n-1) in 
Equation 2.12. Likewise, in treating the large sample, g 
can be substituted for s and y for Hs where y and g are 
the mean and standard error of the entire population. 

It is of interest to determine the probability of 
cccurrence of a particular observation. The normal cr 
GauSSian distribution eguation relates the residual cf a 


particular random variable with the probability of its 


Ze 


cccurrence, and is given by 


1 (5%) 
P(v) = e (2.13) 
o V2n 





The plot cf this equation yields the normal distrikuticn 
curve (Fig. 2.4). Tke height cf the curve above the 
vertical axis 1S propcertional to the probability of a 
particular error occurring. 

The probability of a residual falling between any 
two residuals v and v can be computed by integrating 
pquattion 2.13 ae ; 


; ae 
— 
P(y) = " an (2. 14) 
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Figure 2.4 The Normal Distributioa 


This integral is difficult to evaluate analytically 
so tables have been ccmpiled to aid in computations. Fer 
Ms = +o and a = -o, it can be shown that P(v) = 0.6827. In 
cther words, the protability that a particular observation 


will fall within + lo of the mean is 68.27 percent. 


Ze 


Returning to the example of the vessel moored tc the 
offshcre platform, the mean and the standapganerrer™ fcr tue 
observations are easily computed. With this informaticn and 
Equation 2.18%, the precbability of a Gangemeprore sia aig 
within specified limits can be computed. Conversely, [ry 
fixing a probability, the assoclated ameessoreeane Cane 
error can ke computed . In statistical terms, a particular 
observation will-fall within specified limits with a certain 
confidence. 

Actual values of one-dimensional standard errors for 
hydrographic positioring equipment are a subject of dekate 
tetweer manufacturers and users. Some manufacturers of 
microwave positioning eguipment claim standard errors of +1 
meter. Cn the other hand, Munson [Ref. 9, p. 6] states that 
Microwave systems demonstrate accuracies of 3 meters at 
short ranges but show larger errors at ranges of 15 km and 
greater. NOS assumes a 3-meter standard error in all of its 
Short-range accuracy ccmputaticns. It is apparent that 
further study is needed to adeguately define the nature of 
€Lrors associated with electronic positioning equirment. 

Waltz [Ref. 13] performed an extensive study to 
determine the pointing error of a Wild T-2 theodolite. His 
results showed that the pointing error associated with tnis 
instrument under hydrcgraphic survey conditions was akout 


1.3 meters and was independent of distance. 


2. Two-Dimensional Errors 


The intent of this paper is to apply statistical 
methods developed by cthers tc a hydrographic data set 
containing two-dimensional errecrs which are defined by two 
Tandcm variables. Lengthly and complex derivations are not 
presented. Burt [Ref. 3] and Heinzen [Ref. 2] show adequate 
derivaticns of formulas associated with two-dimensional 


errors and can be referenced for full details. 
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The followirg assumptions are made concerning two- 


dimensicnal errors associated with intersecting LOP's: 


i. The random errers of each LOP are normally 
distri buted. 
ii. Systematic or bias errors have been removed fron 
the observaticns. 
iii. The intersecting LOP's are coplanar. 


1v. The error LOF's are parallel to the exact LOF's. 


In develcping a usarle matheratical model for accuracy 
determinations, the fcur assumptions hold to a high degree 
for all hydrografhic fesitioning geometries. 

Consider again the vessel noored to a fixed crfshcre 
platform. Assume two ranges are measured from two different 
shore control stations at the same time and that the Lrarge 
readings are uncorrelated. The observation orf this pair of 
ranges ig repeated many times. After a large number of 
observations, the means and standard errors of the indi- 
vidual ranges are determined. Suppose the mean rangés, or 
the actual LOP's, intersect at an angle of 909 and that the 
computed standard errcers are equal (g =o )- If each data’ 
ae (iY, ) 1s plotted, the spread of points about the mean 
GeordiNdtes Tesults in a circular cluster (Fig. 2.5). A 
meget CeNSity Of points occurs near the intersection of the 
mean ranges and the density of points decreases outward fron 
the intersection of the mean ranges. 

In this special case, which is called a circular 
Homm@at Gistribution, the probability of a point falling 
within a gepecified radius, R, from the intersection or the 


mean rances 1s 


(omnes 


Sl 


where oe = = = and ws defined aS thenectrewtapeseomeere 
error. Using Equaticn 2.15, R can begcenmputedea eta ag 
P(R), cr conversly, F(R) can be computed by fixing R. 
Letting R = Oo = o. = omme then P(R) = 0.3935. In other 
words, 39.35 percent cf all errors in a circular normal 
distribution are not expected to exceed the circular stan- 


darpd e€rrer [Ret. 12, Eps 25-207. 
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Figure 2.5 Circular Normal Distribution 
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In the case where the two uncorrelated LCP's inter- 


sect at an angle other than 90° or ae ae Peune GCONtC Wins sere 


a 


egual density are ellipses centered about the point defined 


Ey the irtersecting ICP's (Fig. 2.6). The two-dimensional 


probakilaity density function becomes [Ref. 1, p. 136] 
K2 
1 ia (2. 16) 
Eee = ta,c, ° 
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Two Uncorrelated LCP's 


Figure 2.6 Error Ellipse Formed by 


where 
ae is the residual in the direction of the semi-major axis 
of the error ellifse, 


Vy is the residual in the direction of the semi-M1inor 


axe, 
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Oo is the standard error in the direction of the semi- 


major axis, 


o.6C6Uuis) the standard error in the direction of the seni- 


Miner axis, 


and 
v2 oy? 
x 
K2 2 + 2.17 
oe Oy (eZ ) 


The scluticn of Equation 2.16 with values o£ K Eor different 
P's yields the results in Table I [Ref. 12, p. 23]- For a 
39.35 percent probability, the axes of the ellipse are 
1.0000q and 1.0000 go ; for a 50 percent probability, the 
axes Ae 1.1774 O, i 1.1778 oe = 


TABLE I 
Values of the Constant K 





PROBABILITY K 
39435% 1.0000 
50.00% Ve lia 
63.21% 1.4142 
90.00% 2-2-1460 
99.00% | 3.0349 
99.78% | 3.5900 


— 


| 
| 


The error ellipse can te used for accuracy computa- 
tions Ey developing relationships for g ando in terms cf 
the initial informaticn 9 , o , and 55. pocdueen {Ref. 10, 
p- 1213] gives the P gaaloeaets equations for independent I10P's 
relating these quantities: 


Oo 


2. l 2 + 92 + (To? + a2) - 4sin’B a70% 4 
x  2sin7B UCR a Sona ai Baia } (2.18) 
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and 


ee - (Tar Fa)? = F5In7B a7 
oy @sin7B { a pics 4sin°g a*oa } 2 righ 


In these equations, g is assumed to be tne acute angle 
between the LOP's. 

In certain special cases, the above equations take 
Cn more ganageable fcrms. In range-range and aziauth- 
azimuth positioning it is often assumed that on = or = qd. 
Fauations 2.18 and 2.19 then reduce to 


PM 2 ; 
x 2sin(s8) 7520)) 


and 


pice V2. ‘ j 
y 2cos (348) (2.21) 


In the Concentric range-azimuth case, gc #* go , ard ®B 
1 2 
eguals 909. Eguations 2.18 and 2.19 then simplify to 


oC . 
eae (aes 


and 
¥, 2 (22a) 


SUES c > om and o > ss 

The case for ccrrelated LOP's is more complex. fThke 
calculaticn of 9, and oe involves a coordinate transforma- 
tion from a linear skewed coordinate system to an uncorre- 
lated rectargular cocrdinate system. The following 
discussicn is taken from Heinzen [Ref. 2, pp. 49-53]. 

Assume a hydregraphic position is established by the 


mMiGemGecCETOneOr two Correlated LOP*s (rig. 2.Ja). LO? 1 and 
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Figure 2.7 Coordinate Transformations for Correlated ILOP*s 


LOP zZ are the coordinate axes in the skewed coordinate 
systen, with standard errors oF and — The semi-major and 
semi-tiror axes of tke error ellipse are not coincident with 
the skewed coordinate system axes. The correlation cceffi- 
cient hetween the two LOP*s is Pe Assume co > O° 

The standard errors and correlation coefficient in a 
correlated rectangular coordinate system with axes A and B 
must now be determined. A coordinate transformation from 
the skewed system to the correlated rectangular system trust 
Fe made yielding the standard errors along the new coordi- 


nate axes (Fig. 2.7b) 


oo, 2 2 2 (2.24) 
of = =z (of + 2 Oo 0 COSB +0 - 0 
a  sin’sg ( i ey 1 2 ») 2 


St 


and 
come (ees) 


mae Ccrrelation coefficient in the correlated rectangular 
system is 


0 oO 0 2 = ~bs 
=(—cos8 +p ){1l1+p (—)cosB + (—*)cos*g} 
oS 2 12 a _ 


Pab (2.26) 


To determine ae and oat a second coordinate transformation 

must be performed fren the correlated rectangular system to 
an uncorrelated rectangular system with axes X and Y (Fig. 

ze/C). The semi-major and semi-minor axes of the error 


ellipse are then 





(0) a aa 
o, got laa aas 
and 
2 2 2 
C= Oe >t One 
y eae (Z. 28) 
When a = 0, these equations Fecome identical to the 


Simplified versions 1n Bowditch [Ref. 10}. 

The orientaticn of the semi-major and semi-mincr 
axes relative to the intersecting LOP's is the third fparan- 
eter which fixes the error ellipse. The angle 8 (Figs. 2.6 
and 2.7) is measured counter-clockwise from LOP 1 to the 
semi-majcer axis of the error ellipse [Ref. 11] and is given 
Ly 

o7sin(28) + 20 a ao sin(8) 
Laer 22 


@ = i arctan { St 


: o7cos(28) + 20 «o o cos(8) + o (2.29) 
1 12ely 2 2 
For the special case cf es os and p = 0, 
2 
a 5 (are) 
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The orientation of the error ellipse Anvangenenogonal ceor, 
dinate system can be represented by adding or subtracting ¢g 
to the orientation of LOP 1. Care must be taken on deter- 
Mining the quadrant cf the outcome. As a general rule, the 
error ellipse always lies within the acute angles fcrmed Ly 
the intersecting LOP's. 

The orientaticn and dimensions of the error ellipse 
provide a useful index for evaluating the accuracy of a 
hydrographic position. Its greatest attribute is that it 
accurately represents the error distribution about the 
intersection of two ICP's in terms of a fixed probability. 
It is interesting to examine the variation in the relative 
dimensicns and orientations of error ellipses as they vary 
in a Lange-range configuration with 96 =o = 04 Flg< Ze cj 
The dimersions of the ellipses are Spee ct Dye GUatilecrs 
2-20 and 2.21 and omlandes abe funnel lons woes one 
fix edea. There Ronen the ae eens oz the ellipses remain 
constant alorg a contaur Of constant @ ; only tne 
orientation changes. A line of constant 8 is a circle 
which includes staticns L and fF. Note that the dimensicns 
of the ellipses for 8's of 30° and 1509 are identical. the 
ellipses about the 90° angle of intersection contour are 
circles and represent the strongest possible positions in 
this scheme. With varying 8's, the directional nature cf 


the distribution can te noted. 


3. Circular Precision Indexes 


Although the error ellipse gives a true representa- 
tion of the error distribution about a hydrographic fesi- 
tion, its use has certain drawbacks. The characteristics of 
the ellifse must be specified ky the three quantities og , 
oe and @. A Single figure for evaluating the pesnenoaet 
accuracy cannot be used. Greenwalt [Ref. 12, p. 26] states 


that when g and g are not equal, a circular error 
x y 
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Around a Range-Range Systen 


pses 


Ele 


ELITrocr 


distributicn can be substituted for the ells —~ticabedisti mae 
tion. This substituticn can be Satiside toa) =: ogee uc. 
analysis within certain oe ratios. However, when this 
ratio is small the distcrtion introduced by the circular 


distrikuticn may beccre nisleading. 
a. Root Mean Square EIror 


The terms radial error, root nean Square Error, 
and das are identical in meaning when applied to two- 
dimensicnal errors [ Ref. 10, p. 1229Jj. The tern das is 
defined as the square root of the sum of the squares cf tte 
standard errors along the major and minor axes of the errcr 


ellipse. That is 


(Zea) 
where O.. ang GO. are Given Dy EQUations 2. lowandez. ice 
A more direct form of 2.31 is Given by [Reise ce | 


= 1 2 Z 
dans > sing V°, 7% (2-22) 


for uncorrelated LOP‘'s. For range-range and azimuth-azinctth 


positioning, with ae = a =o , Equation 2.32 reduces tec 
d° « V2 0 
rms sing (2.33) 


For range-azimuth positioning, 8 = 90° and Equation 2.32 
kecones 
a (2.34) 


The mere general forg of Equation 2.32 for both correlated 


and uncorrelated LOP's [Ref. 2, p. 59] is 


l a 
rms sind \o2+ ce 29 9 3 cos8 ( ) 


where pP.,is the correlation coefficient. 
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Ah error circle with a radius of one d can be 
constructed about the intersecting LOP's (Fig. 2.9). Two 
d ns is the radius cf the error circle obtained using two 
times the values of o aha oo Ie EQUatTORn 2e3t. FOL an 
elliptical error distribution, the probability associated 
Wet d Specific value of a. varies as a functicn of the 
eccentricity of the error ellipse (Table II). The rgrek- 
ability associated with one ae varies from 63.2 percent to 
68.3 fercent, while the probability associated with two dg. 


varies between 95.4 percent and 98.2 percent. 
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Pigure 2.9 The d Error Circle 


NOS uses dj, aS an accuracy Specrrication. 
Umbach {Ref. 14, p.j 4-25] states that super high frequency 
direct distance measuring systems would be used only when 


Pace va vchot ds 1s less than or equal to: 
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is 0.5 mm at the scale of the Survey for scales of 
1:20,000 and smaller, 
ii. 1.0 mm at the scale of the survey for 1:10,000 
scaie surveys, Or 
iii. 1.5 mm at the scale of the survey for scales of 
1:5,000 and eidrger- 

The major advantage of using G ng as a precision 
index is its ease of computaticn. Some hydrographers draw 
analcgy between the varying probability associated with one 
dist (€2.2 percent to 68.3 percent) and the fixed prer- 
ability associated with a one-dimensional standard errcr 
(63.3 percent). En Lact, ae has very little statistical 


meaning. The obvious problem with using ¢d aS a precision 


rms 
index is the varying frobability associated with the error 
circle. Fer this reason Greenwalt [Ref. 12, p. 31] reccn- 


mends against its use. 


| 
| 
3 


TABLE il 
Probabilities Associated With dt, 


| 
| 
PROBABILITY | 
og oO LENGTH CF | 
vf x 1 Ss 1 dans 2 oone i 
0.0 | ele | 12000 02600 0.954 
Oe lie O To005 Uecoz 0S 5 | 
Or2 | 1.0 | Teo | OO 82 wel Oe oon | 
is oO ol 1.042 i QO.676 {02961 { 
0.5 1.0 1.118 0.662 { 0.969 | 
Ob T20mal 1. 166 | Os650 | 0.975 
Quy 120 -o dees 26 i. O56 4 te 1 20Res7 | 
Omg 1.0 leeza 02635 02250 i 
Oe 10 | 1.345 Os | e991 | 
10 Ae 1.414 | 0 S652 0.982 ! 


is 


L 
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b. Circles cf Equivalent Probability 


Burt [Ref. 3] presents a method for translating 
ellipses of equivalent probability into circles of eguiva- 
dent frokability. Tec utilize this method, it is first 
necessary to compute the eccentricity of the error ellifrse, 
c, by the equation 

0 
ca 
O. UZe410) 


where C4 Sys 

Harter [Ref. 15] compiled Tables III and IV 
which are taken from Bowditch [Ref. 10, p. 1215]. Harter's 
data are given in terms of the eccentricity, c, a parameter, 
K, and a probability, P. The parameter, K, when multipiied 
Dy. os gives the value of the radius, R, of the circle of the 


corresponding probability shown in Table III. That is, 
R= Ko, (237) 


The probatility of a point falling inside a 
circle cf specified radius can be computed by entering Tarle 
III with c and K as arguments. Given a fixed probability, K 
is determined by entering Table iV using c and P as argu- 
ments. The radius of the prokability circle is then 
computed using Equaticn 2.3/7. . 

Using confidence ellipses has certain advantages 
cver confidence circles of equal probability. First, the 
directicnal nature of the true error distribution is act 
represented in the ccnfidence circle method even though beth 
methods give an accurate measure of confidence. Seccnd, the 
area of the confidence ellipse is always less than or equal 


to the area of the ccnfidence circle. The area of a 
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TABLE iV 


Radii of Circles Given c and P 


0.8 | 





| 
7300 | 1. 15035 | 


8 
© 





0. 67449 0. 68199 0. 70585 0. 74993 0. 80785 0. 87042 0. 93365 | 0. 99621 | 1. 05769 1. 11807 | 1. 17741 
1. 15473 1. 16825 1. 19246 1. 23100 1, 28534 | 1. 35143 1.42471 ! 1.50231 ; 1. 58271 1. 66511 
-9000 1. 64485 1. 64791 1, 65731 1. 67383 1, 69918 5; 1. 73708. 5 1. 79152 1. 86253 | 1. 94761 | 2. 04236 2, 14597 
-9500 1. 95996 l. 96253 1. 97041 1. 98420 | 2. 00514 | 2. 03586 | 2.08130 | 2. 14598 | 2. 23029 | 2. 33180 | 2. 447705 
97350 2. 24140 2 24365 2. 25053 aez0200 2. 28073 | 2. 30707 | 2. 34581 2. 40356 | 2. 48494 2.58999 | 2.71620 
.9900 2. 57583 2. 57778 2. 58377 | 2. 59421 ! 2. 60995 | 2.63257 | 2. 66533 | 2.71515 | 2.79069 | 2.39743 | 3. 03485 
= ' { 
-9950 2. 80703 | 2. 80883 2. 81432 | 2. 83289 | 2. 83830 2. $5894 | 2. £8859 2. 93347 3. 00431 | 3.11073 : 3.25525 
9975 3. 02334 2. 02500 3.03010 | 3.03898 3.05234 | 3.07144 | 3. 09871 3. 13969 3. 20556 3.31099 | 3. 46164 
.9990 3 29053 3. 29206 | 3. 29673 3. 30489 | 3. 31715 | 3. 33464 | 3. 35949 3. 39647 | 3. 45698 | 3. 55939 | 3. 71692 
confidence ellirse is 
2 
A= Koo 1 
e Y (2.0) 


where K is the appropriate protability conversion factor 


(Table I). The area of the $0 percent confidence circle is 
Serene (2.39) 


woes R 25 given by EGuation 2.37. For a condition whére 
oP = s = 3 meters, and 8 = 309, the area of the 90 fercent 
confidence ellipse is 261 sguare meters, while the area of 
the ccnfidence circle is 587 square meters. For both stan- 
dard erxrcers equaling 10 meters and 8 = 30°, the 90 percent 
confidence ellipse has an area of 921 Square meters and tke 
confidence circle has an area of 2894 square meters. Fron 
an oferaticnal perspective, the difference in areas kretween 
ellipses and circles have significant implications which 
will te discussed in Chapter V. 

The following examples are presented to demcn- 
strate methods fcr ccmputing the parameters cf error 
ellipses ard confidence circles for several hydrografhic 


positicning geometries. 
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Example J 

A vessel is conducting a hydrographic survey 
using range-range geometry. Tne two .tOr es egencrarcd 2) 
microwave transmitters have standard errors of o = s: neters 
and o = 4 meters. The angle cf intersection 8 ae the 
yeseedl is 309. assume the LOP's are uncorrelated. Ccmpute 
the probability that the vessel"s position Wali ece weer 
circle of 10-meter radius with the center at the intersec- 
tion cf the LOP's. 

Recalling Eguations 2.18 and 2.19, the values of 
0, and o are found tc be 9.79 meters and 6.14 meters, 


respectively. From Eguation 2.36 
Oo 


c= — = 0.633 
Oo 
x 
and frem Egquation 2.37, with R = 10 meters, 
K = 1.032 


Entering Tatle III and using interpolated values for c and 
K, the probability that the vessel's position will be witkin 
a circle of 10-meter radius centered at the intersecticn of 
the LCP's is . 

P = 53.2% 


£xampwe 


In» 


A vessel is conducting a hydroyraphic survey 
uSing range-azimuth cecmetry. The range LOP generated Ly 
the micrcwave transmitter has a standard error of 3 meters. 
The aZimuth LOP determined by theodolite observation has a 
standard error of 1.3 meters at all ranges. Compute the 
radius of the 90 percent confidence circle at the vessel's 


BOs? toms 


In the range-azimuth case 8 = 90° and the LICP's 
are unccrrelated. Tkerefore, 
6 = go = 3.0 meters 
1 x 


and 
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0 = o = 1.3 meters 
Zz 
Soe y 
Oo 
Coe ee oS 
0 
x 
Sablemi¥Y 1S entered with the values of P = 0.9 and c = 


0.432. The value for K is found to be 
= Weed] 
Uomeeeogdatlionezss/, thesradivs of the 90 percent prob- 
apiiaty circle is.fourd to be 
R = 5.14 meters 
Tne probability that the vessel's position will be within a 
circle of 5.14-meter radius centered at the intersection of 


Ene MaGe"s 15.90 —perc ent. 


Example 3 

A vessel is conducting a hydrographic survey 
Mei GgeayeerlLolic-hyperbolic gecmetry. The hyperbolic LC? 
generated by the 1.6-MHZ electronic positioning system hés a 
Standari error of 0.05-lane on the base line. The correla- 
PeonmmcOcerlcClent (5 ) between the two LOP*s is known to tLe 
O.4. Compute the radius of the 30 percent confidence circle 
Meche vessel's position. 

The rectangular plane coordinates of the naster 


(4), two slaves (G and PF), and the vessel's position (P) are 


X COORDINATE Y COORDINATE 
(m) (m) 
R ieee Sat eso 8 
G 208,679.1 08,540.4 
M Zh 7332 Dios att 
P Daa 25 169, 264. 2 


4] 


Given the frequency cf 1.6 MHz, A’ = 187.37 meters fren 
Equation 2.3. The lare width aleng the base line is a = we 
= 93.€8 meters from Equation 2.5. Using the law of ccsines 
from plare geometry, the subtended angles oe and a are 
32.479 and 43.259, respectively. The angle of intersection 
of the two hyperbolas at P is 37.86° from Equaticn 2.7. The 
lane widths at P are a 254.19 meters and w = 335.06 
meters from Equation 2.6. The standard errors of the green 
ic) and red (o ) hyperbolas, respectively are o = w = 


2 1 eg base 
16.7 meters and oe Ww O = 12.7 meters. These standard 


errors are ina 1G eaeweered coordinate system and must Le 
transfcrmed to an unccrrelated rectangular system. fFren 
Eguations 2.18 and 2.19, the values of om and O, are 36%¢ 
meters and 1zZ.7 meters, respectively. The correlaticn coeéf- 
ficiert in the correlated rectangular system (0-4) is then 
0.737 £rcem Equation 2.26. The semi-major and semi-minor 
axes in the uncorrelated rectangular system are 38.1 meters 
and 8.3 meters, respectively, from Equations 2.27 and 2Z.2&. 


The eccentricity is 
0 
c= 2 = 0.218 
x 


Table IV 1s entered with the values of P = 0.9 and c = 
0.218. The value for K is found to be 
K = 1.6602 

From Equation 2.37, the radius of the 90 percent proktakility 
Circle is found atoesc 

Rk = 63.3 meters 
The probability that the vessel's position will be within a 
circle of 63.3-meter radius centered at the intersecticn cf 
the LCP's is 90 percert. 


We 


Titi. EXPERIBENT DESIGN AND IMPLEMENTATION 
The coals of this chapter are to demonstrate that hydro- 
graphic fositioning accuracy can be classified based or tke 
radii cf 90 percent ccnfidence circles determined by using 
Eurt's method and to show that, based on the same criteria, 
accuracy predictions can be made for survey planning 


purposes. 


As. CATA ACQUISITION EROCEDURES 


The data used for analysis and prediction consisted of 
Tange-range, azimuth-azimuth and range-azimuth survey infor- 
mation. The data were acquired by Naval Postgraduate Schcol 
(NPS) students ina Hydrographic Sciences ccurse. Although 
the ccurse was structured as a training exercise, the data 
acguisiticn procedures utilized were nearly identical tc 
those which are fracticed by NOS. 

A total of 453 hydrographic positions were recorded 
during the survey of a nearshore area in southern Monterey 
Pave Calitcepia. Of the positions used for analysis, 292 
were range-range, 81 were range-azimuth, and 80 were 
azimuth-azimuth. All survey information was recorded Ly 
hand in scunding volumes. The vessel used was a 36-foot 
Uniflite with a fiberglass hull and twin engines. The 
survey was conducted cn October 28, November 16, 23, and 30, 
1983. Electronic cecntrol and calibraticn stations used for 
the survey included USE MON 1978, MUSSEL 1932, BEACH LAB 
WIG2Z PUG neehEY ANGRICAN CAN COMPANY STACK 1932, MONTEREY 
RADIG #5 HAT GON KMBY MAST 1962, MONTEREY HARBOR LIGHT 6 1978, 
and MCNTEREY BLUE LIGHTHOUSE (Fig. 3.1). With the excertion 
of MCNIEFEY BLUE LIGHIHOUSE, which is a low-order fosition, 
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Figure 3.1 Hydrographic Survey Area 


all stations are of third-order or better and are published 
in the National Geodetic Survey Data Base. 

FOr azimuth-azimuth and range-azimuth positioning, 
aZimuths were measured with a Wild T-2 theodolite. On 
November 16, range-azinuth information was acguired Ly 
locating the theodolite over station MUSSEL and initialing 
ch USE MCN. The initial direction was eheeked by smenting 
on KMBY MAST. AZimuth-azimuth positions were acguired cn 
Novemter 23. A theodciite was set over USE MON and an 
initial direction was to MUSSEL. A second theodolite was 


set at MUSSEL using USE MON for the inttiaisdirectaon. 
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Range information was recorded using a Racal Decca 
Irispender system, a microwave system commonly used for 
nearshore, line-of-sight survey work. On October 28 and 
Novemker 30, range-range data were recorded Ly setting 
remote units over stations BEACH LAB and MUSSEL. Before and 
after the survey, the ranging system was calibrated cver the 
fixed Fase line USE MCN to MUSSEL. Daily checks in the 
survey area were made to determine if the system was working 
properly. This was accomplished by maneuvering the survey 
vessel tc a point where two known navigational ranges irter- 
sected. One navigaticnal range was formed by stations 
MONTEREY AMERICAN CAN COMPANY STACK and MONTEREY RADIO 
STATICN KMEY MAST. A second navigational range was formed 
Ey staticns MONTEREY EARBOR LIGHT 6 and MONTEREY BLUE 
LIGHTHOUSE. 

Irack control for range-azimnuth and range-range [csi- 
tions was accomplished by steering the vessel alcng range 
arcs. The spacing between range arcs for most lines was 
planned to be 40 meters. Distance between positions alcng a 
sounding line averaged approximately 200 meters. The 
azimuth-azimuth lines were controlled by steering a magnetic 
compass heading. 

The data acquired under training conditions contained 
several deficiencies that would normally not be tolerated. 
For example, the quality of the line steering was generally 
Foor; the vessel wandered off the arc more than 10 meters in 
several instances. The guality of the sounding lines run 
uSing azimuth-azimuth control was extremely deficient; the 
position plot of these lines show a jagged path by the 
vessel. Under normal hydrographic procedures, these pcesi- 
tions wculd be rejected. Since the intent of this study is 
to demonstrate accuracy analysis technicues, these deficien- 
cies prove to be inccnseguential; the acyuired data are 


adequate to demonstrate the concepts. 
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NALYSIS 





IV. RESULTS AND DATA 


baenlA PROCESS LNG 


Autcmated processing of the positional survey data was 
done on the NPS IBM 270/3033AP computer system. Grafhic 
displays were constructed using the Display Integrated 
Software System and Plotting Language (DISSPLA) developed by 
the Integrated Software Systems Corporation (ISSCO) 

[Ref. 16]. All computer programs involved in data 
processing were written in the WATFIV programming language. 

Ccmputations were made in an X-Y coordinate system Lased 
on a Modified Transverse Mercator (MTH) projection. A MTM 
projecticn is essentially the same as a Universal Transverse 
Mercator (UIM) projection, the only difference being that in 
a MIM prejection a central meridian is picked near the 
survey area instead cf being fixed at a particular meridian 
[Ref. 727]. 

The central meridian, controlling latitude, and false 
€asting values define the cocrdinate system used fcr compu- 
taticns. The central meridian for the projection was chosen 
to be longitude 1219 52" 30" W which is approximately the 
mean lcngitude of the survey area. The controlling lati- 
tude, tke distance ir meters from the equator to a reference 
latitude, was chosen to be 4,050,000 meters. A false 
€asting cf 5,000 meters was chosen as the value of the 
X-coordinate at the central meridian. 

Three shore contrcl stations were used in the acgquisi- 
tion of survey data. The geodetic positions of these 
staticns were converted to the X-Y coordinate system (Table 
V) using program UCOMFPS, which is a hydrographic utility 
package available to students at NPS. 
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TABLE VY 
Coordinates of Control Stations 





STATICN NAME GECDETIC CCORD. MTM COORD. 
MON 369 36" 04.685" N | Y = 1982.43 om. 

1710 52% 35.900" wW | X = 4853.36 r. 

MUSSEL 369 37" 18.151" N { Y = 4247.42 mn. 
1719 54" 11.628" Wy X = 2474.75 on. 

EEACH LAB 32€9 36" 05.571" N | Y = 2009.86 om. 
tz 168 52°02. doe | xX = 9914.75. 0. 
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Be. ACCURACY ANALYSIS OF HYDROGRAPHIC POSITIONING DATA 


Tke cbhjective of this section is to illustrate how the 


accuracy of hydrografhic positioning data can be classified 


using Burt's method of circles cf equivalent proprability. 


The radius of the 90 rercent confidence circle was ccmurputed 


for €ach pesition; it frovides a quantitative measure cf 


repeatarle accuracy. 


For sutsequent accuracy computations, the follcwing 


assumpticns were made: 


sles 


asia 


aa 


LVve 


The standard error for the microwave ranging 
system used in the range-range and range- 
azimnuth computations is 3 meters. 

For azimuth-azimuth and range-azimuth positions, 
the pointing error of the theodolite is 1.3 
meters at all ranges. 

The two LOP's involved in all types 

cf positioning are independent ( ao 


The data are free of systematic errors. 


Raw range and azimuth data were hand logged into a data 


file for processing. A modification of program UCOMPS was 


2 


used to compute X-Y ccordinates of all positions.  Basea cn 
cecmetric relaticnshigs discussed ecariter earucies of Titer. 
secticn cf the LCP's vere then computed for rangée-range and 
azimuth-azinuth points. The angles of intersection for all 
rangée~-azimuth positicrs are 909%. 

The range-range and azimuth-azimuth data were then 
passed tc WATFIV subroutine PRCB (Appendix A). As ingut 
parameters, the subrcutine accefts two standard errors cf 
the LCP's and the corresponding angle of intersection. Tke 
output parameters include the semi-major and semi-minor axes 
cf the 90 percent confidence ellipse, the radius of the 930 
percent confidence circle, and the areas covered by both 
figures. 

Subrcutine PROB uses a linear approximation to determine 
the value of the function K for varying values of the eccen- 
tricity, Cc, in Burt's method. A linearwinterpolaciren sya] 
performed by first taking the eleven discrete values of c 
and K for a probability of 90 percent from Tabie IV and tken 
constructing a series of relationshirs for K as a function 
cf£ c (Takkle VI). 

Values of the radii of 90 percent confidence circles for 
rangeé-range data were plotted at their respective fpositions 
(Fig. 4.1). The arcs of circles Connecting themtwercomscL 
staticns BEACH LAB and MUSSEL represent lines of constant 
intersection angle (209°). Of the ranye-range data set, 
position 848 (Appendix B)--coordinates X = 4119.01, Y = 
4735.C7--was found to have the smallest radius (strongest 
Fosition) of 6.4 meters and an angle of intersection cf 
90.29. Position 137--coordinates X = 3345.86, Y = 
3873.34--represents the weakest position with radius value 
of 15.3 meters and an angle of intersection of 26€.7°. 

The positional accuracy degrades rapidly as the inter- 
secticn angle approaches 309; the 30° arc represents a line 


of constant 13.7 meter radius. Within 400 meters of the i300 
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Z 
| | 
| TABLE VI 
| Linear Approxigations for K as a Function of c | 
| Tnterval of ¢ Linear eee | 
Function for kK | 
| 0.0 - 0.1 | K = .0306c + 1.64485 | 
l 0.1 - 0.2 l K = .0940c + 1.63851 
ee he =921652c + 1562427 ’ 
| 0.3 - 0.4 ioe > seer n 250773 
| 0.4 - 0.5 | K = .3790c + 1.54758 | 
05 = 0.6 | K = .5444c + 1.464832 | 
0.6 - 0.7 K = .7101c + 1.36546 | 

pee 08 K = .8508c + 1.26697 
0.8 - 0.9 K = .9475c + 1.18961 | 
| 0.9 = 1.0 | K =1.0361c + 1.10987 | 


intersection arc, the radius varies between 8 and 15 méters. 
The radii values charge slowly in the vicinity of the 
minimus value of 6.4 meters which corresponds to an angle of 
intersection of 90°. 

The radii of 90 fercent confidence circles asscciated 
with the azimuth-azinuth positicns acquired using control 
staticns USE MON and YUSSEL were also plotted at their 
respective fositions (Fig. 4.2). The standard errers cf the 
LOP's are assumed to re 1.3 meters; the resulting imprcved 
accuracy is evident. The maxinpum value of the 90 percent 
confidence circle radii is 8.7 meters at position 
Pew =—coordindtes X = 4327.25, ¥ = 2878.39--which corresfonds 
to an angle of intersection of 159.8° (or in terms of the 
supplement, 20.29). Fosition 682--coordinates X = 4611.20, 
Y = 4421.29--represents the strongest position recorded 
during the survey with a 90 percent confidence circle radius 
of 2.8 meters and an angle of intersection of 91.09. 

Again, the rapid degradaticn of accuracy is noted 
approaching 8 = 150°. The arc of the 150° intersecticr 


angle represents a corstant radius of 5.9 meters. Discrete 
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values along the are corfirm this qualitatively. A large 
area of strong positi¢cnal accuracy SUpeOw G3 geieeitca sa ee 
8B = 90°. Numerous values of 2.8 meters are present near 
the tcp cf the plot. 

Using the assumpticns stated at the beginning of this 
section, the values fcr all radii of 90 percent confidence 
circles for rangé~azinuth positions are 5.1 meters. This 
computation was carried out in Example 2 of Chapter II. 
Since this case is trivial, the data are not displayed 
graphically. 

Positioning data were also classified based on tke 
parameters of the 90 fercent confidence ellipse. WATFIV 
program ELLIP (Appendix C) was used to generate the farame- 
ters cf the 90 percent confidence ellipse for range-range, 
aZimuth-azimuth, and range-azZiputh positioning data. The 
fFrogram waS initialized by entering the coordinates of the 
control stations and standard errors of the LOP's. The f£1x 
number, hydrografhic position coordinates, and angle of 
intersection were then read in from a data file. Sukroutine 
PROB was called to ccmpute values for Ko an Gu haga. 

The angle of orientation of the major axis of the 
ellirse, measured cleckwise from north, was then computed. 
For range-range and azimuth-azinuth positions, the LCE 
generated from the left control station was used as the base 
LOP. For range-azimuth positions, the LOP formed by the 
theodolite was used as the base LOP. First, the orientation 
of the base LOP in tte coordinate system was determined. 

The orientation of tke major axis of the error ellipse rela- 
tive to the base LOP (6) was then computed using Eguaticn 
2-29. By adding or subtracting 9 to the orientaticn of the 
base LOP, the orientation of the major axis of the error 
ellipse in the coordinate system was determined. This angle 
takes on values from 09 to 1809. Appendix D consists cf the 


confidence ellipse classification scheme for range-range, 


azimutk-azimuth and range-azimuth data. Forty positiors for 
€ach positioning gecmetry are listed for comparison to the 
classification scheme presented in Appendix B. 

Appendix B lists the data by position number, X-Y¥ ccor- 
dinate, angle of intersecticn, and radius of the 90 percent 
confidence circle. Appendix D lists the data by positicn 
humber, X-Y coordinate, angle cf intersection, Ko , Ko , and 
angle of orientation for the 90 percent Ponicenice aes 
These afprendices are similar to hydrographic survey data 
bases and demonstrate accuracy classification schemes based 


on the two criteria. 


C. ACCURACY PRELTICTICNS 


bhemeverail positional accuracy of a survey can ke 
contrclled by computing accuracy values before data acguisi- 
tion is kegun. For example, if the hydrographer is using 
radii cf 90 percent ccnfidence circles aS an accuracy 
critericn, the minimunr allowable angle of intersecticn for 
two LCP's can be computed for meeting specifications. the 
Pyeute sOfmenc survey areca Way aliow the flexibility te 
change system geometry to maximize accuracy dt a specific 
location or to maximize the area covered with a given accu- 
racy. By waking acctLracy computations before acquiring 
data, the hydrographer may also have the option of deciding 
what tyre cf positioring system is to be used to meet accu- 
racy reguirements. 

The ccnstruction cf reliability contours is one method 
to display the expected positional accuracy. Reliability 
contours, lines cf constant repeatable accuracy which are 
functions of the system geometry and standard errors cf tke 
rositicrirg equipment, can be constructed about shore 
staticns using the radii of 90 percent confidence circles 


critericn or the less desirable ay value. 


a 


Consider the equations that have been developed in 
Chapter II for the determinaticn of radii of 90 percent 
confidence circles using Burt's method. Fcr uncorreiated 
LOP's in a range-rance or azimuthn-aZimuth systen, the 
repeatable accuracy cf a hydrographic pcsition is a function 
only cf the angle of intersection, assuming the standard 
errors of the LOP's are constant throughout the survey area. 
The Iccus of points which define a constant angle of inter- 
secticn for two LOP's ip a Tanje-trange or aZi1 aie a aa 
system is a circle which passes through both control 
staticns. Given the coordinates of the two control 
stations, the equations of these circles can be determined. 

Construction of reliability contours involves several 
Simple trigonometric relationships (Fig. 4.3). Let IR fe 
the line connecting the two shore control stations L and & 
in a rance-range system. The length of line LR is b. the 
circle through both stations defines a line of constant 
intersection angle for two LOP's. The radius of the circle 
is cr. The distance € 1S measured along the perpendicular 
rkisector of the line IR to the center of the circle at 
Fornt O(h,k) and is given by | 


= b 
e " tang (4.1) 


Knowing e and the radius r, the coordinates of point O can 


re ccmputed. The equation of the circle is then 
r2 = (x - h)? + (y - k)* (4.2) 


These two equaticns were used to generate reliability 
contcurs for display on a computer graphics terminal. Using 
Eurt's method, the angles of intersection of two LOF's were 
computed for discrete values of radii of 90 percent confi- 


dence circles. Reliability ccntours about stations EFEACH 
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LAB and MUSSEL for a range-range systen ocr ma = 3 meters) 


were constructed (Fig. 4.4). Using rau 4 9, X= (ccore 
dinates were generated for points laying on different reli- 
ability circles. A curve-fitting subroutine in the LISSPIA 
library was used to generate the circles through the 
computed peints. The 13-meter accuracy contour corresponds 
to an angle of intersection of 31.6°, while the 7-meter 
accuracy contour corresponds to an angle of intersection of 
67.99. The best achievable accuracy of the system is 6.4 
meters at 909. 

For comparison purposes, reliability contours were 
constructed about BEACH LAB and MUSSEL for azinuth-aZinuth 
geometry (o, zh = 1.3 meters). The increased accuracy of 


this ccenfiguration is evident (Fig. 4.5). The 3-meter 
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contour corresponds tc an angle of intersection of 69.4¢ 
while the 6-meter contour corresponds to an angle cf inter- 
section of 29.69. The best achievable accuracy at an inter- 
secticn angle of 90° is 2.8 meters. 

A second scheme was used tc display accuracy predicticns 
for the twe positionirg methods. Given the coordinates cf 
EFEACH LAB and MUSSEL, a series of discrete points spaced 800 
meters apart, were generated throughout the survey area. 

The values for the radii of 90 percent confidence circies 
were then computed at each point with the use of sukroutine 
FROB. Figures 4.6 and 4.7 illustrate this predicticn 
scheme. These figures present the same information as 
Figures 4.4 and 4.5 in a different manner. The 30° angle of 


intersecticn contour is shown on both figures. 
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Figure 4.4 


GENERATED RELIABILITY CONTOURS 
AZIMUTH—AZIMUTH: BEACH LAB—MUSSEL 
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Figure 4.5 Reliability Contours: AZimuth-Azinuth Geometry 
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ACCURACY PREDICTIONS AZ—AZ 
RADIL OF 90% PROBABILITY CIRCLES 
STATIONS BEACH LAB AND MUSSEL 
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Figure 4.7 Azimuth-Azimuth Point Accuracy Predicticn 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. ACCURACY SPECIFICATIONS 


interpretation of the 1982 IHO positioning standards in 
terms cf 90 percent ccnfidence circles yields some inter- 
esting results with respect to present day survey fractices. 
For example, for a 1:10,000-scale hydrographic survey, NOS 
usually uses microwave fpositicning systems in a range-range 
mode, and assumes a standard error of 3 meters for each LCP. 
Surveys are frequently conducted between the 30° tc 150° 
angle of intersectior limits. Using the 90 percent ccnfi- 
dence circle critericr, the radius of the circle should act 
exceed 10 meters. Hcewever, the radius value for £8 = 30° and 
150° is 13.7 meters. The values of KO and Ko for the 30 
percent confidence ellipse are 17.6 and 4.7 neters, resrec- 
tively. To meet the SO percent critericn for a 1:10,000- 
scale survey, the 8 limits should be 42° to 138°, 

Azimuth-azimnuth fesitioning is accurate enough for 
1:5,00C-scale surveys, using B limits of 35° to 1459, 
assuing a standard error of 1.3 meters for each LCP. With 
the standard error assumptions used for range-azimuth, the 
90 percent radius is 5.1 meters for all positions. Given 
the uncertainties of the standard error figures, it is 
rational to assume that range-azimuth positions can reet tne 
S-meter accuracy standard for 1:5,000-scale surveys. In 
fact, lrange-azimuth fositional accuracy can exceed azinmuth- 
azimuth accuracy wher the later's 8 is less than 359. Fora 
3-meter 0 range-range configuration, it is impossible tc 
meet 1:5,000 specifications with any 8. 

As a general guideline, the 30° to 150° angle of inter- 


secticnh limit is a gccd rule to use for uncorrelated LOF's. 
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However, as mentioned for 1:10,000-scale surveys in a range- 
range mode (o= 3 meters), this rule does not always hold. 
On the other hand, it is possible to have g's of less than 
309 and still meet specificaticns. For example, azinuth- 
azimuth fositioning can theoretically be used for 8's of 180 
to 1629 for a 1:10,00C-scale survey. However, the eccen- 
tricity cf£ the error ¢€llipse is so small that the distcrtion 
introduced by using ccnfidence circles can become 
misleading. In view cf this, eccentricities of less that 
0.2 should not be used. 

Using the 90 percent radius criterion, a table has teen 
assemkled iilustrating the g limit for various positioning 
geometries at different survey scales, uSing assumed stan- 
dard errers (Table VII). The information in Table VII 
illustrates that the 209° to 1509 8 limit need not be fixed. 
The 8 limits should vary based on the scale of the survey 


and the frecision of the positioning eguipment. Accuracy 


-1 meters for Ee = 3 ando 1238 
2 


; TABLE VII | 

| 8 Limits for Surveys | 

| Survey 90% R-R R-R AZ-Az | 

Scale Radius (o = 3) o = 10) (o> gles) { 

B Limit Linit ee Lasmant 

| (m) (deg) (deg) (deg) | 
132,200 | Ze | -=s==>- 2... __- liom ss 

Tee | Delo | maaan | SaaS Je 14S | 
121G7008 1050 42-132 (ee 23> 15 7% 

i 13 26,0007 (5220.0 |..21-153) Ae ip 23-1572 | 

1:40,000 | 40.0 | 23-siSa rae) 35-145 | -23qalo7 = | 

% ee Tami t For Oe | 
Note: 90% radii cf all range-azimuth positions are 

| assumed to be | 

{ 

4 
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figures as a function of 8 for uncorrelated LCP's have keen 


compiled using standard errors of 
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azimu eh 


(Table VIII), 3.0 meters for range-range short-range 


(Table IX), and 10 meters for range-range medium-range 
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DeeeUeko POR ACCURACY FIGGRES 


NCS is currently developing the Shipboard Data System 
Betas tithe, a hydrcgobaphic data acquisition and 
processing system which will replace the present HYDROLCG,/ 
HYDRCFLOI system. SLE IIL will revolutionize data acguisi- 
tion and processing technigues with the capability tc 
perforg high-speed calculaticns and display color graphics. 
With this increased ccmputer potential, data manipulaticns-- 
such aS accuracy comrutations--can be performed. 

Each position in a survey can be given a guality figure 
based on the radius of the 90 percent confidence circle. 
This figure is sufficient for non-critical positions of 
ordinary hydrographic data. Critical positions are thcse 
which are determined for significant features (i1.e., wrecks, 
least depths, rocks, and other potential hazards). fFcr 
these fositions, the farameters of the 90 percent error 
ellipse can be computed, as well as the radius of the 90 
percent’ confidence circle. 

Many schemes can te envisioned for the use of an accu- 
racy figure. For example, suppose the position of a 
submerced pile was determined Ly range-azimuth geometry in a 
prior survey. The radius of the 90 percent confidence 
Gurere sic theheo.! meters (Ex. 2, Ch. Ii). The chartirg 
agency now wishes to relocate the pile to determine if it 
Still exists and is still a hazard to navigation. In lcw 
water visibility, a ccmmon technigue used to resolve such an 
item would be to send divers down over the reported fosition 
and ccnduct a circle search. One diver remains at the 
reported position, hclding a line, while the other diver 
swims a circumference holding the other end of the line. 
Theoretically, if the line is about 5 meters long and a hang 
does ret occur, it is 90 percent certain that the file has 


keen remcved. For a higher confidence, the iine is 


a4 


lengthened. In an investigation such as this, it is advi- 
Sable tc be conservative and use the maximum length of line 
which is operationally feasible to provide coverage cf an 
area as large as possible. The radius of the 90 percent 
confidence circle gives the hydrographer a rough figure fcr 
answering the questicn: Does the submerged pile exist? 

Knowing the farameters of the error ellipse couid fe 
useful fcr conducting wire-drag, wire-sweep, and side scar 
sonar oferations. Fcr a position obtained with low freci- 
Sion positioning eguirfment, the search to relocate a 
submerced feature could cover a large area. Knowing the 
parameters cf the error eilipse could reduce tke area, tite, 
and effort of the search. The search pattern could te 
planned to cover the desired confidence ellifse. 

With the quantification of accuracy, a decision must fe 
made concerning how mruch confidence is needed to delete a 
certain feature from the chart after a search has Leen made. 
The 90 fercent confidence level may be too low, whereas tte 
$5 or 99 percent level may suffice. A balance must Le main- 
tained retween confidence of disproval and time and effcrt 
spent on the search. 

Accuracy predicticns in the form of reliability contours 
can ke displayed using computer graphic terminals. These 
displays will contrikute to the efficient planning of 
surveys to meet specifications. Given the survey area, tte 
availarkle ccntrol, the positioning methods, and the preci- 
sion of the positionirg equipment, the hydroyrapher can plan 
the accuracy of the survey before it is conducted. The 
Survey area and the available control may be such that there 
is flexibility to change control stations to optimize accu- 
racy over an area of critical importance. This informaticn 
can be displayed graphically and plans for the survey can be 
made acccrdingly. Likewise, given an accuracy limit, such 
aS a 10-meter radius cf the 90 percent confidence circle, 


the area tc be covered at that accuracy can be wnaximized. 
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Many variables exist when considering accuracy reguire- 
ments fcr a hydrographic survey. In general, higher accu- 
racy means more time, money, and effort. Azimuth-azimuth 
geometry is the most accurate method of positioning analyzed 
in this thesis. This method involves at least two people 
asnore and good ship-to-shore communications. Currently, 
NOS acquires these data manually, which minimizes the speed 
that the vessel can cferate and adds to processing tire. On 
the otker hand, a survey using a medium-range system needs 
little shcre support and the data acquisition is autcmated. 
Accuracy predictions kelp keep a balance between accuracy 
and effort. If the desired accuracy is attainable using a 
range-range system instead of an aZimuth-azimuth systen, 
then tke chcice is okvious. 

Hydrcgraphic pcsitioning in the future will be dcmirated 
ry two methcds. For cffshore surveys, the Global 
Positioning System (GFS) is expected to give positional 
accuracy to 10 meters or better. GPS is a satellite posi- 
tioning system currently being deployed by the Department of 
Defense and will provide near worldwide coverage for users. 
Since tke full constellation of 18 satellites will not te 
operational until 1988, it 1s not yet known if the expected 
accuracy of 10 meters will be met. Nearshore surveys may 
use multiple LOP's for establishing hydrographic pcsiticns. 
The principle of least Squares 1S applied to redundant 
cbservations yielding the most probable position. Fer toth 
GPS and least Squares positioning, confidence ellipses and 
circles can be determined, although the technigues invcived 
are much mere complicated than those presented in this 
thesis. 

The accuracy classification scheme presented in this 
thesis is predicated cn the elimination of systematic 
errors. Much work is needed in identifying the sources of 
systematic errors asscciated with hydrograrhic positioning 


equipmert. 
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PARAMETERS 


SUERCUTINE FOR 90 FERCENT CONFIDENCE CIRCLE 
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{C.GT.0.8)) THEN 
GT.0.7)) THEN 
-GT.0.6)) THEN 


INEAR INTERPOLATION 


—.0.9)) THEN 


TY OF ELLITES: 
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PLS been 


U 


C 


=0.7101*C +1. 36546 
ELSE IF ((C,1E, 0,6) -AND. (C- 
=0.5444*C+1.46488 
ELSE IF ( C-LE.0-5) AND: (C. 
=0.3790*C+1.54758 
ELSE IF ((C-1E20.4) . AND, (Cs 
=0.2535*C+ 1.5977 
ELSE IF ((C.1E.0.3) AND, (C. 
=0.1652*C4#1.62427 
ELSE IF ((C-TE. 0.2) -AND: (C. 
=.094*C#1. 6385 
ELSE 
B=0.0306*C+#+1. 64485 
END IF 
RADIUS=B¥*SICX 
SGX90=2.146¥*S1IGX 
SGY90=2. 146 #SIGY 
CIRAR=3.1415626*RADIUS**2 
ELAR=3. 14 15926*SGX90*SG190 
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RANGE-RANGE ACCURACIES 
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RANGE-~RANGE ACCURACZES 
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APPENDIX C 


PEOGEAM FOR 90 PERCENT CONFIDENCE ELLIPSE PARAMETERS 


: PROGEFAE NAME: ELLIE 
C DESCRIFTION: COMPUTES ORIENTATION, 90% SIGMA-X, 90% SIGMA-Y, 
c FCR ERROR ELLIPSE ABOUT A HYLDECGRAPHIC 
c POSITICN ESTABLISHED BY RANGE-RANGE, AZIMUTH- 
: AZIMUTE OR RANGE-AZIMUTH POSITION 
Seo obueCns NECHOLAS E. PERUGINI 
C LT. NOAA 
c NAVAL FOSTGKRADUATE SCHOOL 
: DATES SEPTEMBER , 1984 
e 
f Pirorerd REAL ~*~ (A-H,C=-Z) 
Cee Nit PATI ZESVALUES: FOR RANGCE-=RANGE AND AZIAUTH-AZIMUTH: 
C “KL AND YL ARE COORDINATES OF LEFT STATION 
c Sun oND TReAnRE COORDINATES OF REGHT STATICN 
c ~SIGi AND SIGR ARE RESPECTIVE STANDARD ERRORS 
: ASSCCIATED WITH EACH LOP 
C FOR KRANGE-AZIMUTH: 
c =XipoNE Yr ARE CCORDINATES OF OCCUPIED STATION 
c —loeetowo LOMA Or TREODOLITE Lop 
c -SiGR Is SIGMA OF RANGE LOP 
CK OK A OR eK RK A RK KK Ae eK a KK A kK eK KO KKK RK Ke 
XL=4914.75 
YL=2009. 86 
: SIGL=3.0 
XR=2474.75 
YR=4247.42 
P SIGR= 3.0 
ect erecenwe GAd PON CCEPFICITENT: USUALLY ZEKG FOR R-E&, 
c R-AZ, AND AZ-AZ 
- PT=3.141593 
C ENTER INDICATOR TO TELL WHAT ATM? OF CATaAwtS ENTERING PROGRAY 
C IND = 1 KANGE-RANGE 
C IND = 2 AZIMUTH-AZIMUTH 
C IND = 3 RANGE-AZIMUTH 
CK A OM Me Me Me Ke eK KK KK eK KK OK ek ke KOK OK 
IND=1 
oie BEE Eh Auth cu eA AAR p Praha ee ee ee 
C INC IS A TOGGLE WHICH CHECKS FOR BETA GREATER THAN 90 DEG. 
C NOTHING IN PROGRAM SHOULD BE CHANGED FROM HERE ON 
C 
Conta ternevAtA FROUW CATA FILE: IFIA,= Fix NUMBER 
ec EX =e COCOnRDUNALTE Ore nYDRO@ POSITION 
e Eee =r eooORDENATE, OR RYDR@. POSITION 
ey ID = ANGLE OF INTERSECTION IN DEGREES 
i Sotho se rte = SI9, TELLS PROGRAM TO SIOP REALING 


10 CCNTINUE 
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rsd 
= >< 
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F7.2, 5X,F8.4,3X,F5-2) 
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WOW © 
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Oo ~~ Jd 
O * 
hor 
* &) 
Ow 
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Inne 
He 8 by bdO 


oO 


HOU 


DEGREES: TOGGLE TRATED Cn Ic Cie 


tH 
*mhtex<ae: 
+e» 


Otte ¢ 


OQ 
On © 
Fb 

tr 


CHANG 


SiS FCR COME GraATZons 
L BE FIZED Sie RESP BC Ree Oe ee ee 


> hrf fet 


RTH BETWEEN HYDRO POSITION AND 
LL DOr ace PEIWEEN O= MeOeey tr 


Hi 
GE AZIMUTH DETERMINATION. 
0 40 
) 


Lea 
GREES 


ANAAANAA AA 


THEN 
PI-ATAN ((PY-YL) /(PX-XL) ) 
ATAN ((PY-YL) /(XL-PX)) 


ELSE 
THEN 

= ATAN((YL-PY)/(PX-XL)) 

PI-ATAN((YL-PY) /(XL-PX)) 


b> 
ES 
md 
ttt 
to 

lI 


GC TO 60 
gene FIXING FOR AZIMUTH-AZIMUTH POSITIONS 


4 CCNTINUE 


C 
c 
C 
ash = ATANS ((PX-XL)/(PY-YL)) 
ALPHA = PI-ATAN((XL-PX)/(PY-YL)) 
IF (PX.GE.XL) THE 
ALPHA PI-ATAN ( (PX- XL) / (YL-PY) ) 
ALPHA = ATAN((XL-PX) /(YL-PY)) 


AZIMUTH EQUALS THETA FOR RANGE AZIMUTH CASE, ASSUMING 
TEEODCLITE SIGMA IS LESS THAN RANGE SIGMA 


IF (IND.BQS3) GoMBiog7o 


MPUTING THETA, THAT IS THE ANGLE OF ROTATICN FRCM 


ANOAA QANAAN 


BEGIN CO 
LEFT LCP 
60 CCNTINUE 

B1=SIGL** 2¥*SIN (2*BETA) +2*RO* SIGL*SIGR*SIN (BETA) 
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B2=SIGL** 2*COS (2*®BETA) +2*RO*SIGL*SIGR*COS (BETA) +SIGR¥*2 
IF (ABS (B2).L1.0.0001) E2=.0001 
B2=E1/B 
C COMFUTE ROTATION ANCLE FROM LEFT LoP 
TH=0.5*A TAN (B3) 
; 90 CCNTINUE 
€ 
C DEFINE SEMI-MAJOR AXES ORIENTATION IN TERMS OF 0-180 LEGREES 
c ROTATICN, CLOCKWISE FROM NORTH 
C RANGE-RANGE CASE 
IF(IND.EQ.1) ‘TEEN 
IF INC- EQ. 1) THETA=ALPHA+TH 
IF (INC. EO. 0) THETA=ALPHA-TH 
: END I 
C AZIMUTE-AZIMUTH CASE 
(IND.EQ.2) TEEN 
IF INC. EQ.) THETA=ALPHA+tTH 
IF ({INC.EO. 1) THETA=ALPHA-TH 
END I 
c 
C FIX RCIATION ANGLE FROM 0-180 DGREES 
70 CCNIINUE 
C CCNEITICN FOR RANGE-AZIMUTH DATA 
IF (IND. EQ.3) THETA=ALPHA 
IF THETA LT~O;) THETA= EI+THETA 
TE(THETA.GT.PI) THETA= THETA-PI 
C DEG IS THE SEMI-MAJCK ELLIPSE AXIS ORIENTATION IN DEGREES 
DEG=57. 295779¥* THETA 
C COMEUTE 90% SIGMAX AND SIGMAY OF ERROR ELLIPSE 
CATT. PROB (SIGL SIGR,RO, ID SGX90,SGY90 RADLUS, ELAE, CIRAR) 
WRITE (7 1 0) SFIX,PxX,PY,TD 5Gx50 SG150, 
100 FCRMA {ix-1 BK FT 3x 67 P22 708° 5 Pied. Seer 4 oa. 
* 3X, FS. ) 
6C’T0 1 
900 CCNTINUE 
STOE 
END 
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